The Aharonov-Bohm effect in graphene rings 
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This is a review of electronic quantum interference in mesoscopic ring structures based on 
graphene, with a focus on the interplay between the Aharonov-Bohm effect and the peculiar elec- 
tronic and transport properties of this material. We first present an overview on recent developments 
of this topic, both from the experimental as well as the theoretical side. We then review our recent 
work on signatures of two prominent graphene-specific features in the Aharonov-Bohm conductance 
oscillations, namely Klein tunneling and specular Andreev reflection. We close with an assessment 
of experimental and theoretical development in the field and highlight open questions as well as 
potential directions of the developments in future work. 
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I. INTRODUCTION 

The Aharonov-Bohm effect^- is a fundamental phe- 
nomenon of quantum interference related to the trans- 
mission of particles through a closed loop pierced by 
a magnetic flux. Besides it's fundamental significance 
for quantum theory, it's importance for applications in 
mesoscopic interferometric devices such as the electron 
Sagnac gyroscope* is omnipresent. The effect was origi- 
nally observed in metal rings in 1985^ and later also in 
carbon nanotubes£ In graphene, an atomically thin two- 
dimensional carbon allotrope which was first isolated in 
2004;£ the Aharonov-Bohm effect is expected to exhibit 
unusual behavior due to the peculiar electronic proper- 
ties of this material, in which charge carriers at low en- 
ergy behave effectively as massless Dirac fermions, giving 
rise to a number of (pseudo-)relativistic effects such as a 
Berry's phase ir— & (for reviews on graphene see Refsliol- 
T2I ). In the following, we review recent developments on 
the Aharonov-Bohm effect in graphene nanostructures, 
first from the experimental side and, later on, we elabo- 
rate more on theoretical aspects. 

Experimental progress. The first experimental re- 
alization of a graphene ring structure was reported 
in 2008^ In this work, the authors investigate the 
Aharonov-Bohm oscillations in diffusive single-layer 
graphene as a function of the magnetic field which is 
applied perpendicular to the graphene plane in a two- 
terminal setup. They find clear magnetoconductance 
oscillations with the expected period corresponding to 
one magnetic flux quantum $0 = h/e on top of a 
low-frequency background signal due to universal con- 
ductance fluctuations which are present in any disor- 
dered, phase-coherent mesoscopic device. Increasing the 
temperature T gives rise to thermal averaging of the 
Aharonov-Bohm oscillations, and the authors find that 
the oscillation amplitudes decay as T -1 / 2 , as commonly 
observed in metal rings»« 

The authors further observe two unusual features in 
the recorded data. First, they find indications of a lin- 



ear relationship between the oscillation amplitude and 
the overall ring conductance. Such a behavior has nei- 
ther been observed in metal rings, nor in semiconductor 
heterostructures, with the exception of ReffTol where a 
similar effect is seen. The authors speculate that tunnel 
barriers which may be present in their device could be re- 
sponsible for the observed behavior; however, a detailed 
theoretical analysis has yet to be done. 

A second peculiar feature is the significant increase of 
the oscillation amplitude at strong magnetic fields close 
to the onset of the quantum Hall regime. This increase is 
strong enough to make the second harmonic — i. e. oscil- 
lations of period $o/2 = h/2e — visible in the frequency 
spectrum. Such a behavior was also observed by an- 
other group in subsequent experiments with smaller rings 
and higher visibility in a two-terminal as well as a four- 
terminal geometry and was attributed there to scattering 
on magnetic impurities^ — an explanation derived from 
corresponding observations in metallic rings that is not 
compatible with the observations made in Refill, where, 
the authors instead speculate that the increase of the 
oscillation amplitude may be due to orbital effects orig- 
inating from a potential asymmetry in the arms of the 
ring; however, this assumption could not be confirmed in 
subsequent numerical calculations ji£ 

In Refill, the authors also introduce additional tun- 
ability into the graphene ring device by applying a side 
gate potential to one of the ring arms. In subsequent 
experiments^ the same group systematically investi- 
gates the influence of such side gates in a four-terminal 
geometry (see Fig. Q]) in the diffusive regime and find 
phase shifts of the Aharonov-Bohm oscillations as a func- 
tion of the gate voltage as well as phase jumps of tt at 
zero magnetic field — direct consequences of the electro- 
static Aharonov-Bohm effect (which is more feasible in 
graphene than in metal rings due to the low screening 
of this material) as well as the generalized Onsager re- 
lations. The authors of Reffl6l further speculate on the 
presence of edge disorder indicated by the fact that the 
various charge carrier trajectories around the ring, which 
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FIG. 1: (a) Scanning force micrograph and (b) magnetoresis- 
tance of the four-terminal (Sl/2, Dl/2) graphene ring struc- 
ture with two side gates (SGI/2) investigated in Reffla (Fig- 
ures adapted from Reffla). 

can be derived from the frequency distribution of the 
magnetooscillations, do not cover the full area of the ring 
arms. 

Further experiments on graphene ring structures in- 
clude the local oxidation nanolithography using atomic 
force microscopy 19 and antidot arrays on epitaxial 
graphene films^; in the latter setup, universal conduc- 
tance fluctuations are suppressed since the sample size 
exceeds the phase coherence length, while Aharonov- 
Bohm oscillations are still visible due to the small size of 
the antidots. As a graphene-specific feature, the authors 
also observe increased visibility of weak localization due 
to intervalley scattering on antidot edges. Shubnikov-de 
Haas and Aharonov-Bohm effects on thin graphite sin- 
gle crystals with columnar defects were investigated in 
RefsEm in a four-probe measurement, and a signifi- 
cant contribution of surface Dirac fermions ("graphene 
on graphite") as well as evidence supporting the the- 
oretical prediction of edge states was found. Another 
multi-terminal measurement of Shubnikov-de Haas oscil- 
lations in monolayer graphene relating the Landau level 
separation between electrons and holes with the trans- 
port gap in the density of states is given in Refl23l: the 
poor visibility of the Aharonov-Bohm oscillations therein 
is attributed to a low phase coherence length. 

Theoretical progress. On the theoretical side, there 
is a great variety of topics which have been addressed, 
including the valley degree of freedom characteristic for 
graphene, particular device geometries and edge symme- 
tries, resonant behavior and transistor applications, as 
well as the role of interactions. Common quantities ex- 
pressing the influence of these aspects include electronic 
properties such as spectrum and persistent current, as 
well as transport properties such as conductance and 
noise. 

In the pioneering work on the topic of Aharonov-Bohm 
rings made of graphene^ it was shown — both analyti- 
cally using a circular ring in the weak intervalley scatter- 
ing limit in the continuum model, as well as numerically 
using a hexagonal ring structure with zigzag terminated 
edges and strong intervalley scattering — that in the con- 
fined geometry of a graphene ring structure an applied 
magnetic flux gives rise to a lifting of the orbital de- 



generacy in a controllable fashion, which is manifest in 
conductance and persistent current, even in the absence 
of intervalley scattering. The magnitude of the lifting of 
valley degeneracy and it's dependence on details of the 
geometry of the ring have subsequently also been dis- 
cussed in Refill. The effect was further observed in nu- 
merical magnetotransport simulations in the closed (or 
weakly coupled) circular ring geometry described by the 
tight-binding modelii for the case of smooth mass con- 
finement, where intervalley scattering is suppressed. In 
this work, simulations have been done both for ballis- 
tic and diffusive regimes, and up to the quantum Hall 
regime where Aharonov-Bohm oscillations are found to 
be suppressed. 

Perfectly shaped ring devices with clean, well defined 
edges such as the hexagonal ring structure mentioned be- 
fore^ have been addressed in various studies, exploiting 
graphene-specific features of such ideal nanostructures. 
A detailed numerical study of the influence of shape and 
geometry, edge symmetries and corner structures on the 
electronic structure in the presence of a magnetic field 
reveals for example the edge state anticrossing and there- 
fore gap opening due to the coupling of states localized 
at the inner and outer edges of the ring as well as the 
crucial role of corners in zigzag or armchair edge ter- 
minated rings^ For instance, the corners in an ideal 
hexagonal ring with zigzag edge termination in a mag- 
netic field introduce intervalley scattering, as was con- 
sidered in a supercell approach within the tight-binding 
model—; in this work, also a peculiar dependence of the 
spectrum as well as the persistent current on the even 
(resulting in semiconducting behavior) or odd (resulting 
in metallic behavior) character of the number of atoms 
across the ring arm width was found. Similar rings with 
metallic armchair termination have been considered in 
Refs |28l29i ; in such systems, appropriately chosen cor- 
ner junctions exhibit signatures of effective broken time 
reversal symmetry — caused by pseudomagnetic fields — 
at low energies, such as broken particle-hole symmetry 
or a gap in the spectrum that may be closed by the 
application of a real magnetic flux (The generation of 
pseudomagnetic fields in graphene rings under strain (or 
shear stress) has further been discussed in Ref|3(ih. In 
Refl3ll. the authors calculate the spectrum and the per- 
sistent current in ideal diamond shaped graphene ring 
structures either with zigzag or armchair edge termina- 
tion within the tight-binding model; they also encounter 
the even-odd behavior mentioned before, and compare 
their results with previous work on hexagonal rings 31 

Such ideal structures have also been found to be domi- 
nated by resonant behavior in the magnetotransport, for 
instance in ReflU where for small rectangular graphene 
nanorings with perfect edges resonant tunneling through 
quasi-bound states was observed rather than Aharonov- 
Bohm oscillations, which may be tuned by varying geom- 
etry, Fermi energy, or magnetic field. Resonant behavior 
was also observed in RefUTI as well as in in Ref.33j, where 
it was proposed to utilize the electrostatic Aharonov- 



3 



Bohm effect via sidegates — such as already realized ex- 
perimentally in Refll8l — for application in a quantum in- 
terference transistor with high on/off ratio made of a 
hexagonal graphene ring structure with perfect edges, 
where armchair edges were found to be preferable. In 
contrast, zigzag edges in the leads, acting as valley niters, 
and a circular graphene ring structure exhibiting an irreg- 
ular boundary were considered in Ref|34l, where resonant 
behavior was encountered as well; however, the main 
finding of this work was that for opposite valley polariza- 
tion in the leads and appropriately sized rings exhibiting 
higher harmonics in the Fourier spectrum due to multiple 
turns around the ring region, the lowest harmonic is sup- 
pressed while higher harmonics are unaffected. A similar 
geometry, namely a graphene Aharonov-Bohm ring con- 
nected to valley filters which encircles a dislocation, was 
considered in the continuum model in Refl3ol and deco- 
herence properties were discussed. Possible applications 
as quantum interference transistors have also recently 
been discussed in Ref|H for disordered graphene rings, 
where ballistic rectification and negative differential re- 
sistance are observed in the /-^-characteristic; further, 
while for temperatures as large as 150 mK, phonon scat- 
tering is negligible, future work on the effect of electron- 
phonon interaction might be interesting. 

Other aspects of interactions in graphene quantum 
rings have also been addressed in Refl37l in the contin- 
uum model, with a focus on the interplay between valley 
polarization and Coulomb interaction, affecting the val- 
ley degeneracy. This influence was found to be accessible 
through the fractional nature of the periodicity of the 
Aharonov-Bohm oscillations in the persistent current as 
well as changes in the absorption spectrum. An analy- 
sis in the tight-binding model including electron-electron 
interactions further revealed the connection between elec- 
tronic correlations and the spin polarization of the inter- 
acting ground state as function of ring size and number of 
electrons . 38 ' 39 The interaction with the electromagnetic 
held has been discussed in Refill where it was found 
that, in a graphene ring threaded by a magnetic flux, ex- 
citations generated by electromagnetic pulses give rise to 
"valley currents" — in analogy to spin currents. 

The peculiar properties of graphene systems suggest 
to consider the Aharonov-Bohm effect also in more 
graphene-specific ring geometries and transport regimes. 
The electronic properties of monolayer as well as bilayer 
graphene rings in a magnetic field have been addressed in 
Refsl4l|-|43l. either defining the ring geometry by tuning 
the band gap of bilayer graphene or employing a sim- 
plified zero-width ring geometry within the framework 
of the Dirac equation. In the magnetotransport within 
the (Andreev-)Corbino disk geometry in graphene, Lan- 
dau level resonances, the suppression of conductance os- 
cillations away from the charge neutrality point, and 
the crossover to the normal ballistic transport regime at 
large doping and weak fields have been encountered. 44 
Analytical expressions for conductance and Fano fac- 
tor in the magnetotransport of pseudodiffusive graphene 
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FIG. 2: (Color online) Schematic of the graphene ring struc- 
ture (left) and the y-dependence of the on-site gate potential 
V (right) that is applied to the lattice sites on the lower arm 
of the ring. V exhibits either a smooth (black solid line) or 
sharp (red dashed line) profile along y-direction while being 
constant along x-direction. The width of the arms of the ring 
is chosen equal to the width w of the leads. The shaded area 
indicates the region of non- vanishing homogeneous magnetic 
field pointing out of plane. Different edge disorder configura- 
tions are realized by randomly removing sites within the two 
regions of width S between dashed and solid circles (Figure 
adapted from Ref fig!) . 



rings, where transport at the Dirac point is dominated 
by evanescent modes, have been derived in Refl45l 

In the next sections, we will review our previous work 
on two graphene-specific effects in more detail. In Sec. [[[1 
Klein tunneling of Dirac fermions through a potential 
barrier in a graphene ring sy stem similar to the exper- 
imental setup of Refs lfaflg is considered. A related 
effect — Andreev reflection at a graphene-superconductor 
interface in a mesoscopic ring device — is subsequently 
discussed in Sec. IIIII In Sec. IIVI we will conclude and 
provide an outlook on potential future developments on 
the topic of the Aharonov-Bohm effect in graphene. 



II. INTERPLAY OF THE AHARONOV-BOHM 
EFFECT AND KLEIN TUNNELING IN 
GRAPHENE 

In this section, we review the numerical investigation 
of the effect of Klein tunneling 46,47 on the Aharonov- 
Bohm oscillations in a graphene ring on the basis of a 
tight-binding model with nearest-neighbor couplings. In 
order to introduce Klein tunneling into the system, an 
electrostatic potential can be applied to one of the arms 
of the ring (see Fig. [5] for a schematic), such that this 
arm together with the two adjacent leads form either 
a nn'n- or ropn-j unction (n, n': conduction band trans- 
port, p: valence band transport). The former case corre- 
sponds to normal tunneling and the latter case to Klein 
tunneling. Then, the transmission properties strongly 
depend on the smoothness of the pn-interfaces. In par- 
ticular, for sharp junctions the amplitude profile is sym- 
metric around the charge neutrality point in the gated 
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arm, whereas for smooth junctions the Aharonov-Bohm 
oscillations are strongly suppressed in the Klein tunnel- 
ing as compared to the normal tunneling regime. Such 
a setup thus allows for a clear graphene-specific signa- 
ture in Aharonov-Bohm measurements which seems to 
be readily observable. Its physical origin is the quantum 
interference of normal tunneling as well as Klein tunnel- 
ing trajectories through the two arms of the ring. The 
workpresented in this section was originally published in 
Refill 



A. Model 

The calculation is based on the usual tight-binding 
Hamiltonian for graphene 



0*1. 



(i) 



where the second sum runs over nearest-neighbors and 
Vi = V{vi) is a position-dependent on-site potential, tak- 
ing the origin of coordinates at the center of the ring. The 
graphene hopping integral tq ~ 2.7 eV picks up a Peierls 
phase in the presence of a magnetic field, yielding for the 
nearest-neighbor coupling element the expression 



-t exp 



— / A(r)dr 



(2) 



where the line integral is taken along the straight path 
between sites i and j. $o = h/e is the magnetic flux 
quantum, and 



A(r) =-By6(d-\x\)e x 



(3) 



with d = \J R 2 — w 2 /4 is the vector potential giving rise 
to a homogeneous magentic field 



B(r) = V x A(r) = B 6(d - \x\)e z 



(4) 



The system under consideration is a ring-shaped struc- 
ture cut out of a graphene sheet, which is attached to 
two crystalline leads also modeled using the graphene 
lattice structure (see Fig. [2]). Besides the magnetic field, 
the structure is also subject to a gate electrode poten- 
tial V g located on top of the lower arm of the ring. The 
smoothness of the potential interface is controlled via the 
smoothing width w s measured from the lower edges of the 
leads: 

V = i) for y>-w/2, 
V = V g for y < — w/2 — w s , 
< V < V g otherwise. 

In the presented simulations, a cosine-shaped smoothing 
profile is used and < w s < R — 3w/2. 

For a Fermi energy E > 0, together with the adjacent 
leads this lower arm forms either a nn'n- or npn-junction 




FIG. 3: (Color online) Schematic of the influence of the po- 
tential profile introduced by V g on the spectrum of the lower 
arm of the ring. The left hand side shows the normal tun- 
neling case (nn'n-junction) and the right hand side the Klein 
tunneling case (npn-junction) (Figure adapted from Ref l48h . 



for Vg < E and V g > E, respectively (see Fig. [3] for 
a schematic; the Fermi energy is measured relative to 
the charge neutrality point in the leads). Note that the 
setup exhibits a flat potential profile for trajectories along 
the upper ring arm, i. e. a nnn-junction, since there is 
no gate potential applied. This enables a rather large 
transmission through the ring even when the lower ring 
arm is tuned towards the Dirac point, since transport 
through the upper arm always takes place at an energy 
distance E away from the charge neutrality point. 

Transport calculations are done in the Landauer- 
Buttiker formalism for elastic transport at zero tempera- 
ture assuming complete phase-coherence using a variant 
of the recursive Green's function technique^— for the 
coupling of the surface Green's functions of the leads, as 
described in RefUl; the linear conductance of the system 
is expressed through the transmission function Tr (tH), 
where t is the N x N transmission matrix between the 
two leads, i. e. the lower left block of the unitary scatter- 
ing matrix S, 



S 



r a 

t r' 



(5) 



which is obtained from the Green's function of the cou- 
pled system via a Fisher-Lee relation. N is the number 
of propagating modes in the leads. 



B. Results 

In the following, we present transmission properties for 
a ring with R/clq — 300 and w/ao — 60, ao = 0.142 nm 
being the nearest-neighbor distance in graphene, in terms 
of the linear conductance G = 2e 2 /h ■ Tv(vt), where the 
factor 2 accounts for spin degeneracy. Edge disorder is 
applied to the ring by randomly removing sites within 
a width S from the inner and outer edges of the ring, 
respectively (see Fig. We choose S/a = 1.5 in or- 
der to keep the edge of the ring as smooth as possible 
while still allowing for different edge disorder configura- 
tions. Fermi energy E and gate potential V g are chosen 
such that transport always takes place in between the van 
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FIG. 4: (Color online) Magnetoconductance of a ring with 
R/ao = 300, w/ao = 60 at energy E/to = 0.5 and zero gate 
voltage, showing clear Aharonov-Bohm oscillations on top of a 
background due to universal conductance fluctuations (Figure 
adapted from Ref|H). 



Hove singularities located at E = ±r where the density 
of states diverges in the tight binding model of graphene, 
< E < r , < V g < 2 E. In Fig. H we plot the magne- 
toconductance at Fermi energy E/tq = 0.5 and zero gate 
voltage (V g = 0) for a particular ring realization, showing 
pronounced Aharonov-Bohm oscillations on top of a low 
frequency background. The background signal results 
from universal conductance fluctuations (UCF) which are 
typical for phase-coherent mesoscopic devices^ The be- 
havior is in agreement with the observations made in 
Reffl7l. where the authors investigate an even wider mag- 
netic field range up to the quantum Hall regime. In 
Fig. we also show the corresponding frequency spec- 
trum obtained from a Fourier transform of the magneto- 
conductance signal up to B = 1CP 3 Qo/a 2 ,, as well as the 
UCF background signal and the magnetooscillations after 
background removal by means of a high pass frequency 
filter. The contributions to the Aharonov-Bohm oscil- 
lations are centered around (ABa§e//i) _1 ~ 2.3 • 10 5 . 
Using R 2 7r • AB — h/e, this frequency corresponds to a 
mean radius R/ao ~ 270 of interfering electron trajec- 
tories, which perfectly lies within the boundaries of the 
ring. 

By applying a gate voltage V g > to one of the ring 
arms, the magnitude of the Aharonov-Bohm oscillations 
may be modified. A convenient measure of the oscilla- 
tion magnitude is the root mean square (RMS) ampli- 
tude of the signal. Prior to the RMS analysis, the UCF 
background has to be removed from the signal. This 
is achieved by applying a high pass frequency filter to 
the Fourier transform of the magnetoconductance data, 
as indicated in Fig. [5] The retained, unbiased signal is 
squared, and the root of the average over the squared 
signal is defined as the RMS amplitude AGrms- 

In Fig. |5] we show the dependence of the RMS os- 
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FIG. 5: (Color online) (a) Frequency spectrum corresponding 
to Fig. [31 obtained from the Fourier transform G of the magne- 
toconductance G. Besides the low frequency background and 
the fundamental oscillation component, the second harmonic 
is also slightly visible in the spectrum. The dashed line indi- 
cates the frequency limit of the high pass frequency filter used 
for background subtraction (Figure adapted from Ref |48T) . (b) 
UCF background signal (red) and magnetooscillations (blue) 
corresponding to Fig. fJJ after background subtraction. 
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FIG. 6: (Color online) RMS analysis for the setup 
used in Fig. f4j for different smoothing widths w s /ao £ 
{0, 21, 52.5, 105, 210}. Each data point results from an aver- 
age over five realizations of edge disorder. The corresponding 
standard deviations lie between 0.005 -2e 2 /h and 0.015 -2e 2 /h 
but are suppressed for better visibility. For better clarity, the 
spectrum schematics (see Fig. [3| are also included (Figure 
adapted from Ref 14811. 



cillation amplitude AGrms on the gate voltage V g for 
different smoothing widths w s (see Fig. [2]) at energy 
E/to = 0.5, where the average is taken over the range 
< Bciq/Qo < 10~ 3 . Increasing the gate voltage from 
zero towards the neutrality point V g = E not only leads 
to increased potential scattering but also to a reduction in 
the number of accessible propagating states in the lower 
arm of the ring. As can be seen in Fig. [51 the oscilla- 
tion amplitude diminishes and reaches a minimum value 



6 



at the neutrality point. Note that, since the transmis- 
sion through the upper ring arm is not at all affected 
by a gate potential, the overall conductance itself is only 
slightly changed to fluctuate around 2.5-2e 2 //i (see Fig. 7 
in Ref,4&) as compared to values around 3.4 • 2e 2 /h in 
the case of zero gate potential on the lower ring arm (see 
Fig.®. 

For V g < E, the decay of the RMS amplitude towards 
the neutrality point does not depend on the details of the 
gate potential interface. However, in the regime of Klein 
tunneling, V g > E, the oscillation behavior strongly de- 
pends on the smoothness of the gate potential. In case 
of a smooth potential, the partial waves in the lower arm 
have to tunnel through a finite region of low density of 
states, where V ~ E (see Fig. [2J), in order to interfere 
with the partial waves traversing the upper arm. The 
lower arm becomes increasingly penetrable as this region 
gets narrower, until it gets transparent in case of a sharp 
potential. This reflects the usual behavior of Klein tun- 
neling phenomena, where the probability for tunneling 
through a pn-junction depends on the smoothness of the 
pn- interface! 11 ' 47 

The described behavior of the RMS amplitude is ro- 
bust over the whole energy range under consideration, 
except for an increasing uncertainty at lower values for 
the Fermi energy. Although all results are presented for 
zigzag boundary conditions in the leads, the effects are 
independent of a change of orientation of the graphene 
lattice to armchair boundaries in the leads. 

Klein tunneling in graphene thus exhibits clear sig- 
natures in the Aharonov-Bohm oscillations observed in 
mesoscopic rings. In the next section, we will show 
how another graphene-specific effect, namely specular 
Andreev reflection at a graphene-superconductor inter- 
face^ can be identified in such nanostructures. 



III. HOW TO DISTINGUISH SPECULAR FROM 
RETRO ANDREEV REFLECTION IN 
GRAPHENE RINGS 

In this section, we review numerical transport calcu- 
lations of Andreev reflection in a graphene ring system 
threaded by a magnetic flux and attached to one normal 
conducting and one superconducting lead. To this end, 
the Bogoliubov-de Gennes equation for the tight binding 
model using the recursive Green's functions technique is 
solved within the Landauer-Biittiker framework for elas- 
tic transport. By tuning chemical potential and bias volt- 
age, it is possible to switch between regimes where elec- 
tron and hole originate from the same band (retro config- 
uration) or from different bands (specular configuration) 
of the graphene dispersion, respectively. Andreev reflec- 
tion is known to be closely related to the effect of Klein 
tunneling discussed in the previous section^ However, 
different aspects of Klein tunneling have become exper- 
imentally accessible in the last years ; 55 ' 56 whereas spec- 
ular Andreev reflection has not been observed to date, 



although there exist a number of proposals for the ex- 
perimental control 5 - 7 - and detectio n 53 ' 58 ' 59 of this process. 
(For a review on both effects, see RefE.) Here, we re- 
view a novel approach concerning the identification of 
specular Andreev reflection, distinguishing it from con- 
ventional retro reflection, and discuss the advantages over 
previous works in the field. We find that the dominant 
contributions to the Aharonov-Bohm oscillations in the 
subgap transport are of period h/2e in retro configura- 
tion, whereas in specular configuration they are of period 
h/e. This result confirms the predictions obtained from 
a qualitative analysis of interfering scattering paths, and 
since it is robust against disorder and moderate changes 
of the system, it provides a clear signature to distinguish 
both types of Andreev reflection processes in graphene. 
The work presented in this section was originally pub- 
lished in Ref® 



A. Scattering path analysis 

Our approach is based on the observation, that in gen- 
eral, the probability for an incident electron to be re- 
flected as a hole is less than one. This allows for effects 
typical for phase-coherent mesoscopic devices, like uni- 
versal conductance fluctuations or Aharonov-Bohm os- 
cillations 2 in the magnetoconductance. While in nor- 
mal metals, the fundamental period of these oscillations 
is given by the flux quantum $o = h/e, it is half the 
value for Andreev (retro) reflection in conventional met- 
als, due to the charge 2e of a Cooper pair. However, 
this is not true anymore in the case of specular Andreev 
reflection, therefore providing a criterion to distinguish 
between specular and retro reflection. In order to show 
this, we consider the phases due to the magnetic flux 
that are picked up by the various scattering paths. In 
this analysis, we restrict ourselves to the contributions 
up to first order in the sense that we take processes into 
account that involve only a single electron-hole conver- 
sion process, and that contain at most one additional 
round-trip of electron or hole, respectively; higher order 
contributions connected with additional round-trips are 
often times negligible! 13 ' 48 The corresponding paths are 
summarized in Fig. [7J In order to obtain the magneto- 
conductance for the two types of Andreev reflection, we 
sum up the amplitudes as defined in Fig.[7]for the various 
paths coherently: 

i? s ($) Si \s+ + s_ + s' + e i$ + s'_ e" 1 *! 2 

Rr($) = \r+ e i$ + r_ e" 1 * + r' + e 2i * + r'_ e - 2i *| 2 

where s' ± = s' ±e + s' ±h , r' ± = r' ±e + r' ±h , and $ is the 
magnetic flux measured in units of the flux quantum $o- 
Assuming |s| 3> |s'| for any zeroth- and first-order am- 
plitudes, respectively, we obtain 

R s (<P)^R° s + 2Re[(s' + s*+.soS-)e i *} +0[(s')\ (6) 
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FIG. 7: (Color online) (a) Device geometry showing a 
graphene ring structure that is penetrated by a magnetic flux 
<3> measured in units of the flux quantum $o- At the inter- 
face with the superconductor (shaded region), electron- hole 
conversion may occur, (b) The gauge is chosen such that 
each of the eight individual electron (solid lines) and hole 
(dashed lines) paths picks up a phase ±$/2 as indicated, (c) 
Scattering paths for electrons injected from and holes leav- 
ing through the left normal conducting lead; only zeroth and 
first order contributions are included, i. e. terms containing a 
single electron-hole conversion process and at most one addi- 
tional round-trip of the electron or the hole. The paths are 
categorized according to the total phase that is picked up, 
and each path is associated with a corresponding amplitude, 
where first order amplitudes are indicated by a prime (Figure 
adapted from RefmOl). 



where sq = s + + s_ and R® contains contributions that 
are constant with respect to $. Therefore, in the case of 
specular reflection, oscillations of period h/e are domi- 
nant. In contrast, in the case of retro reflection, contri- 
butions of period h/2e are dominant, as expected: 

i? r ($)^i?,° + 2Rc[r + rlc 2i$ ] + 0[r r' , {r'f], (7) 

where again R® contains ^-independent terms and we as- 
sume \r\ ^> \r'\ for any zeroth- and first-order amplitudes, 
respectively. 



B. Numerical model 

In order to test the previous analysis on the basis of 
a microscopic model, we implement the Bogoliubov-de 
Gennes Hamiltonian 61 



H — Ep A 
A Ep — H* 



(8) 



within the tight binding formalism of graphene, Eq. (fl} . 
In this numerical calculation, all higher order contribu- 
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FIG. 8: (Color online) Schematics of the (radially symmet- 
ric) excitation spectrum (lower panel) and surfaces of con- 
stant excitation energy in k-space (upper panel) in the cases 
Ef > e > {retro configuration) and < Ef < e (specular 
configuration). Solid and dashed lines indicate electron- and 
hole-like states, respectively, (hole) states originating from 
the valence band are shaded. The small arrows in the up- 
per panel indicate the direction of propagation of the corre- 
sponding states. Electron-hole excitations are drawn assum- 
ing conservation of k y at the NS interface (Figure adapted 
from RefUJ. 



tions beyond the ones discussed in the previous subsec- 
tion are also taken into account. In Eq. (|8]), we assume 
A,; = A(r,) 6 1 for the superconducting order parameter 
A = J^i \i) (i\ ■ The presence of a magnetic field is cap- 
tured by a Peierls phase in the hopping matrix element, 
Eq. ©. 

The structure of the graphene device under consider- 
ation is schematically shown in Fig. [7] The two semi- 
infinite leads also exhibit the graphene lattice structure; 
superconductivity is induced into the right lead due to 
the proximity effect of a superconducting electrode on 
top of the graphene. We choose to orient the leads to 
exhibit armchair edges and later comment on the reason 
for this particular choice. The whole ring is penetrated 
by a uniform perpendicular magnetic field of strength B, 
described by the vector potential ([3]). The origin of co- 
ordinates is taken at the center of the ring. 

In order to fulfill the mean-field requirement of super- 
conductivity, which demands the superconducting coher- 
ence length £ = Hvf /A to be large compared to the wave- 
length As in the superconducting region^ we introduce 
additional doping into the superconducting region by ap- 
plying a gate potential Vi =V g 8(xi — d). Which type of 
Andreev reflection occurs at the NS interface is then de- 
termined by the excitation energy e (i. e. the eigenvalues 
of Eq. ([8])) and the Fermi energy Ep, as shown in Fig. [8] 



In retro configuration, Ep > e > 0, where Vy ■ v^f 1 < 
for the y-components of the electron and hole velocities, 
both electron and hole traverse the same arm of the ring. 
In specular configuration, < Ep < e, the hole is re- 
flected back through the other arm of the ring, since 
. ,,( e ) > o. In the following, we choose \V g \ ^> Ep, 



u 



8 



justifying the adoption of the step- function model for the 
superconducting order parameter, A, = A0(xj — d)J& 

In order to compare retro (r) and specular (s) configu- 
rations, we will choose e' 1 *' = Ep and = Ep since 
then the states in both configurations exhibit the same 
wavelength and there is the same number of propagat- 
ing modes. We further choose e( r \E]p <C e^ s \Ep^ so 
that for nearly each value of k y , there exist electron- hole 
scattering channels. 

Again, the transport properties of the system are ob- 
tained from the scattering matrix S that is calculated in 
the framework of the Landauer-Biittiker formalism us- 
ing the recursive Green's function technique as in the 
previous section. In the framework of the Bogoliubov- 
de Gennes Hamiltonian, Green's function and scatter- 
ing matrix are parameterized by the eigenvalues e of the 
Hamiltonian ((5J. 

In the following, we will concentrate on the regime 
e < A, in which there are no propagating modes in the 
superconducting lead, so that electrons injected from the 
normal conducting lead are reflected back either as elec- 
tron (e) or hole (h). The scattering matrix thus has the 
structure 



s 



7*ee 'eh 
The r hh 



from which the differential conductance for the Andreev 
processes is given by 

dV~ h l ^ T he r he) 

where the factor 4 accounts for spin degeneracy and the 
quantization of charge in units of 2e. 



C. Results 

In Fig. 121 we show the calculated transmission for a 
ring of width w = 87 ao and outer radius R = 500 ao, 
where ao is the distance between nearest neighbors. The 
transmission function exhibits Aharonov-Bohm oscilla- 
tions on top of a low frequency background which is due 
to universal conductance fluctuations. The position of 
the NS interface is given by d = 400 ao- The chosen 
dimensions of the ring are large enough to exclude finite- 
size effects while still being numerically manageable. For 
the superconducting order parameter, we choose a value 
of A = 0.03 To ~ 80meV, which may appear unrealistic 
at first sight, considered the fact that typical values are 
up to a few meV. However, by making this choice we 
scale the value of the superconducting order parameter 
according to the scale of the system size, such that the di- 
mensionless factor AR/Hvf stays of same order of magni- 
tude, compared with values realized in experiments . 13 ' 18 
Thus, for a realistic system size of R ~ 10~ 6 m, our choice 
of A would correspond to a value of a few meV for the 
superconducting gap. Note that due to these low energy 
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FIG. 9: Differential magnetoconductance for specular (black) 



and retro (gray) configuration for Ep 
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corresponding to 8 modes in the nor- 



mal conducting lead, including all degeneracies (spin, valley, 
electron/hole). The high doping in the superconducting lead 
is chosen such that Ef — V g = 0.5 to in both cases. Other pa- 
rameter values are provided in the main text. The period of 
the dominant oscillation is w 1.8TO _6 ao~ 2 /i/e in specular 
configuration and Bq w 8.8 ■ 10 _7 a ( y 2 ft/e w 0.5 in retro 
configuration. The weak beating pattern in retro configura- 
tion and the asymmetry in specular configuration arise due to 
minor contributions of contrary frequencies (Figure adapted 
from Ref®. 



scales and the rather large spacing of modes resulting 
from the narrow geometry of the electron waveguides in 
such a ring structure, in specular configuration only the 
regime of a low number of modes is accessible. Also note 
that due to strong electron backscattering at the front of 
the hole and at the rough edges of the ring, the average 
value of the differential conductance is much less than a 
conductance quantum, e 2 /h. 

The average radius f of the scattering path is calcu- 
lated according to f 2 irBo = h/ne, where n = 1 (n = 2) 
in specular (retro) configuration and Bo is the (domi- 
nant) period of the oscillation. Evaluating the period of 
the oscillations shown in Fig. El we 

obtain « 420 ao m 
specular configuration and s» 425 ao in retro configu- 
ration. The obtained values lie well within the inner and 
outer radius of the ring and close to the arithmetic mean 
R — w/2 £3 425 ao, therefore confirming the predictions 
obtained from Eqs. ([6]) and ([7]). Minor contributions of 
period h/e in retro configuration and h/2e in specular 
configuration visible in Fig. [9] may arise due to terms 
neglected in Eqs. (0) and (0, scattering off the sharp 
boundaries of the ring structure, and the fact that for 
the electron-hole conversion at the NS interface k y is not 
strictly conserved. 

Another strong evidence that supports our interpreta- 
tion of the two different periods is the breakdown of this 
particular signature that is observed for a shift of the 
position of the NS interface on the scale of the width 
of the ring. Indeed, while in Refill — where a three- 
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FIG. 10: (Color online) Breakdown of the h/e vs. h/2e sig- 
nature for shifted positions of the NS interface, as explained 
in the text. Other parameters and color coding are chosen as 
in Fig. [9] For d = 340 ao (left), in specular configuration one 
observes oscillations of period h/2e as in retro configuration. 
For d = 490 ao (right), contributions of specularly reflected 
holes in retro configuration become important, leading to the 
observation of additional ra/e-oscillations. The value of the su- 
perconducting coherence length is £ = 50 ao (Figure adapted 
from RefH|). 



terminal graphene junction is analyzed — the exact po- 
sition of the NS interface has no effect, it matters here; 
the reason is that £ is comparable or even less than the 
system size, while in Refl62l the superconducting coher- 
ence length greatly exceeds the system dimensions. If the 
interface is too close to the hole region (see Fig. [TU](a) in- 
set), then specularly reflected holes are forced to traverse 
the same arm as the incoming electron. In this case, one 
should observe h/2e oscillations in specular configura- 
tion. In contrast, if the interface is too far from the hole 
(see Fig.[in](b) inset), holes may significantly be reflected 
through the other arm, e. g. due to increased scattering 
at the ring boundaries. This would manifest itself in the 
observation of h/e oscillations in addition to the h/2e 
oscillations in retro configuration. This behavior is con- 
firmed in the observed magnetooscillations, as shown in 

Fig. EH 

Apart from that, the h/e vs. h/2e signature proves to 
be very robust against moderate changes to the length 
and energy scales in the system, such as the extent of the 
magnetic field or the ratio of Fermi wavelength and the 
width of the NS interface. We also tested that the signa- 
ture persists when more propagating modes are present 
in the lead, leading to values of the average conductance 
which are much larger compared to the few-mode situ- 
ation shown in Fig. [5] (see Fig. ITT1) . Additionally, the 
signature is hardly affected by bulk disorder, which is 
a major advantage of our setup. In Fig. [T^l we show 
the magnetoconductance of the system used in Fig. [5] 
with a particular random short-range disorder configura- 
tion, which is realized by applying an uncorrelated, ran- 
dom on-site potential of Gaussian distribution with zero 
mean and width a — 0.01 tq to each site. In addition, 
the NS interface has been smeared out over a distance 
I = 90 ao in this case. The robustness of the effect can be 
explained from the topological nature of the signature: 
since all microscopic scattering paths can be classified 
into just two groups — yielding h/e- or /i/2e-oscillations, 
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FIG. 11: Differential magnetoconductance for d = 430 ao, 
E ( p ] = 0.029 r = e (s) , = 0.0001 r = £ (r) , correspond- 
ing to twice the number of propagating modes, compared to 
Fig. [9] Other parameters and color coding are chosen as in 
Fig. [9] Note that the magnetoconductance signal is strongly 
enhanced compared to Fig. [9] 




B [I0- 4 $ /ag] 

FIG. 12: Magnetoconductance of the system used in Fig. [9] 
with a smooth potential profile (inset) with I — 90 ao and bulk 
disorder of strength a — 0.01 To as explained in the text. The 
h/e vs. h/2e signature still persists. The color coding is the 
same as in Fig. [9] (Figure adapted from Ref |60T) . 



respectively — according to which arm is traversed by the 
quasiparticles, impurity scattering and the resulting de- 
flection of quasiparticles has no adverse effect as long as 
scattering between the groups is weak, while scattering 
within one group may be arbitrarily strong. In addi- 
tion, note that while our description of transport via the 
scattering matrix assumes complete phase coherence, a 
signature that distinguishes retro from specular Andreev 
reflection is assumed to persist also in the case of a finite 
phase coherence length. More specifically, if the phase 
coherence length is on the order of the ring circumfer- 
ence, first-order amplitudes in Eqs. © and (J7J may be 
neglected. Then, retroreflection would still manifest itself 
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in /i/2e-oscillations, while there would be no oscillations 
at all in the case of specular reflection. 

Finally, we add a remark concerning the choice of arm- 
chair boundary conditions in the leads employed in the 
analysis in this section. In a tight binding implementa- 
tion of graphene, there are two simple choices for the ori- 
entation of the leads. Often, zigzag edges are considered 
to represent a generic boundary condition for graphene 
ribbons^ In this case, edge states are present in the sys- 
tem that modify the simple picture provided in Fig. [5] 
by adding additional scattering channels between bulk 
and edge states while removing certain scattering chan- 
nels between bulk states due to the conservation of the 
so-called pseudoparity symmetry that acts like a selec- 
tion rule^ In the realistic limit of metal leads providing 
a large number of propagating bulk modes, this effect 
should be less important. However, for the system geom- 
etry used in the numerical calculations in combination 
with the low energy scales, it may significantly affect the 
observed behavior. In order to avoid this influence, we 
chose armchair boundary conditions in the leads that do 
not provide any edge states. Note in addition, that in 
realistic systems the zigzag-specific effect would also be 
suppressed since the zigzag edge state is not protected 
against disorder when next-nearest neighbor hopping is 
taken into account^ Therefore, we are convinced that 
our results based on armchair edges in the reservoirs de- 
scribe the generic situation for wide leads. 



there is a great variety of aspects covered by theoret- 
ical considerations, and despite the significant amount 
of work done on graphene ring systems, there are still a 
number of open questions drawn from initial experiments 
that remain unanswered so far, as well as a large num- 
ber of theoretical predictions not yet confirmed by cor- 
responding experiments. On the one hand, it would be 
interesting from a theoretical point of view to investigate 
the origin and significance of the seemingly linear rela- 
tionship between conductance and oscillation amplitude 
as well as the significant increase of the oscillation ampli- 
tude at high magnetic field observed in Reffl3l. Further 
analysis on the role of interactions may also be worth- 
while. On the other hand, while some of the theoretical 
models are hardly realizable in experiments at the present 
stage e. g. due to insufficient control over edge properties, 
there are systems, such as presented in the previous two 
sections, that should be experimentally accessible and 
robust, and should therefore allow for the observation of 
graphene-specific features in the Aharonov-Bohm effect. 
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